Weak localization in disordered systems at the ballistic limit 
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The weak localization (WL) contribution to the two-level correlation function, R(uj), is calculated 
for two-dimensional disordered conductors. Our analysis extends to the nondiffusive (ballistic) 
regime, where the elastic mean path is of order of the size of the system. In this regime the 
structure factor, S(t), (the Fourier transform of -R(w)) exhibits a singular behavior consisting of 
dips superimposed on a smooth positive background. The strongest dips appear at periods of the 
periodic orbits of the underlying clean system. Somewhat weaker singularities appear at times which 
are sums of periods of two such orbits. The results elucidate various aspects of the weak localization 
physics of ballistic chaotic systems. 

PACS numbers: 73.20.Fz, 72.15.Rn, 03.65.Sq 
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I. INTRODUCTION 



Interference effects, arising from the interplay of phases 
accumulated along different paths, are particularly inter- 
esting in ballistic chaotic systems. It is due to the hier- 
archy of importance among the classical trajectories in 
these systems: Long trajectories exhibit a universal sta- 
tistical behavior, while short ones constitute the dynam- 
ical fingerprints of the system. The stable (and therefore 
usually also the shorter) the orbit is the stronger is its 
signature. This signature appears both in the wave func- 
tions (the scar phenomenon^) as well as in the statistical 
properties of the energy spectrum of the systemH. The 
purpose of this paper is to study the fingerprints of the 
classical periodic orbits on the nature of interference in 
chaotic systems. 

Our best understanding of quantum interference is in 
disordered systems. In these systems interference may 
lead to localization of the particle in spaceB. If, how- 
ever, the disorder is too weak to localize the particle, 
interference manifest itself as an increase in the return 
probability compared to the classical value. This effect, 
known as weak localization (WL), has been observed by 
measuring the magnetoresistance of metallic filmsu. 

Recent advances in nanostructure technologyB, opened 
the possibility of manufacturing clean mesoscopic sys- 
tems - systems in which the elastic mean free path, I, is 
of order of the size of the system, L. It is natural to ask 
what is the analogue of WL in such ballistic systems? 

Very little is known about this issue, mainly be- 
cause of the failure of periodic orbit theory to provide a 
simple systematic procedure for calculating interference 
(i.e. WL) correctionstl. This failure has been one of the 
main motivations for constructing the. supersymmetric 
nonlinear cr-model of ballistic systemsQ. The hope was 
that this model will produce a WL expansion for bal- 
listic systems analogous to that of disordered systems. 
However, it turned out that WL crucially depends on 
the regularization of the fieild integral, and only specific 
cases could be worked outfl. These are the cases where 
the dynamics is still diffusive or dictated by random ma- 



trix theory (RMT)i. 

Usually one would choose to study the WL signature 
in transport properties, because they are naturally re- 
lated to the experimental data. However, this choice will 
be inappropriate for our purpose for the following reason: 
WL (similar to localization) takes place on a certain man- 
ifold in phase space. For example, in disordered systems 
this manifold is the real space, while in a circular billiard 
with rough boundaries localization occurs in the angular 
momentum spaceEj. In general chaotic systems there is 
no preferred basis, therefore, WL ma;y-appear on a com- 
plicated manifold in the phase spacefill. Yet, transport 
measurements dictate a preferred basis, and may totally 
miss the WL physics we seek to describe. 

Nevertheless, interference effects manifest themselves 
also in the spectral properties of chaotic systems, which 
are independent of the choice of basis. Therefore, in this 
work we shall focus our attention on the WL contribu- 
tion to the simplest nontrivial spectral quantity - the 
two-level correlation function: 



R(uj) = A 2 (p(e + fiw)p(e))-l. 



(1) 



Here p(e) — J2 a ~ e «) ^ s the density of states, A = 
l/(p) is the mean spacing between neighboring energy 
levels, e a , and the averaging, (•••), is over the disorder 
configurations or the energy e. 

To state our problem in this context, consider the den- 
sity of states of quantum system with Hamiltonian hav- 
ing a classical chaotic counterpart. Gutzwiller's trace 
formulallj expresses the density of states, in the semi- 
classical limit, as a sum over the classical periodic orbits 
of the system: 
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where s p (e) is the action of the p-th periodic orbit, and 
A p is the corresponding amplitude depending on the sta- 
bility of the orbit and its periodti 

The traditional way of calculating correlators such as 
(|l|), within periodic orbit ^theory, is to use the so called 
diagonal approximationEja. In this approximation one 
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FIG. 1. (a+b) An illustration of pairs of orbits contributing 
to the diagonal approximation: (a) an orbit with itself, (b) an 
orbit with its time reversed counterpart, (c) The eight-shaped 
periodic orbits associated with the WL contribution to R(u>) 
in diffusive systems. 

replaces a double sum over periodic orbits (such as that 
obtained when substituting (0) in (nT)) by a single sum: 
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where t p = ds p (e)/de is the period of the p-th orbit. The 
rational behind the diagonal approximation is that the 
coherent contributions, at u) = 0, come from pairs of or- 
bits (p,p f ) having precisely the same action. Thus one 
should pair orbits with themselves, p = p' (Fig. 0a), as 
well as with other orbits related by symmetries such as 
time reversal symmetry (Fig. |l|b). In the absence of other 
spatial symmetries, (3 in the above formula is one for sys- 
tems with time reversal symmetry, and two for systems 
which do not have this symmetry. 

The problem of WL in the context of the two-level 
correlation function can be formulated as: How can one 
improve on the diagonal approximation to include inter- 
ference effects systematically? 

In seeking for the solution of this problem it is natu- 
ral to inquire about the situation in disordered systems 
where the systematic interference corrections to the di- 
agonal approximation is the "weak localization" expan- 
sion. The diagrammatic picture of the WL correction 
to R(u>) suggests that the WL contribution is associ- 
ated with pairs of periodic orbits crossing-.themselves at 
some point in space as shown in Fig. ffldiil. Thus along 
one loop the two orbits propagate in the same direction, 
while along the other loop they are in opposite directions. 
However, such orbits exist only in the presence of a non- 
classical scattering potential, and do not have a direct 
analog in the periodic orbit theory. 

Facing this difficulty, in this work, we study WL using 
disorder diagrammatics but far from the diffusive regime, 
i.e. when the elastic mean free path is of order of the size 
of the system. In this case, the disorder is sufficiently 
weak, and traces of the short periodic orbits of the un- 
derlying clean system are still significant. 

We, thus, consider a system consisting of a particle 
confined to move on a two dimensional torus, in the land- 
scape of a random potential, see Fig. |2| The Hamiltonian 




FIG. 2. (a) An illustration of the model used in this paper 
for calculating the WL effects in the ballistic limit. The sys- 
tem consists of a non-interacting electron gas on a torus with 
white noise potential. This potential is sufficiently weak, such 
that the elastic mean free path is of order of the size of the 
system, (b) An equivalent representation of the system as a 
square with periodic boundary conditions. 



of the system is 



(4) 



where p is the momentum of the particle, m is its mass, 
and V(r) is a Gaussian random potential defined by 



(V(r)) = 0, and {V(r)V{r')) 



2ttvt 



(5) 



Here, v = 1/AL 2 is the averaged density of states per 
unit area, and r is the elastic mean free time for scat- 
tering on the potential. This system has been consid- 
ered earlier by Altland and Gefenlla and by Agam and 
FishmanEj, but only in the framework of the diagonal 
approximation. 

In analyzing the results of the above model, it will be 
convenient to use the spectral structure factor defined as 



S(t) 
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(G) 



Using S(t) one can relate the quantum spectral proper- 
ties of the system to the behavior of its classical analog. 
In particular, S(t) form a connection to the periodic or- 
bits of the system: Substituting (|J) in (||) one sees that, 
within the diagonal approximation, the structure factor 
takes the form of a sum over peaks located at times which 
equal to the periods of the classical periodic orbits: 



S(t) 



2hA 



\A p \ 2 6(t-r p ). 



(7) 



It has been noticed by Argaman et al.0 that the right 
hand side of the above equation can be also interpreted 
as |i|p(i), where p(t) is the classical return probability at 
time t. The notion of return probability has been further 
developed by Chalker et al. to obtain a more accurate de- 
scription of the structure factor for diffusive electronstS. 

A disorder potential usually erases the (5-singularities 
of S(t) associated with the classical orbits of the clean 
system. But, if it is sufficiently weak, it will leave traces 
of the them. Indeed, S(t) calculated, in the diagonal apt 
proximation, for weak disorder, shows a series of peaksE£l 
(see inset of Fig. 3). The locations of these peaks 
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FIG. 3. The structure factor of chaotic systems with time 
reversal symmetry. The solid line represents the results of 
random matrix theory. Magnified is the regime where pertur- 
bation theory applies and nonuniversal corrections, which are 
the main focus of this paper, are important. Here we depict 
only the results of the contribution of the "diagonal approxi- 
mation" . The peaks, indicated by pairs of winding numbers, 
are the signatures of the periodic orbits of the clean system 
(see Fig. 4). Both, the Fermi velocity and the system size, L, 
are set to unity. The elastic mean free path, in these units, is 
1/2. tu = 2nh/A is the Heisenberg time. 

along the time axis are precisely the periods of the orbits 
of the clean system. (These orbits are defined by pairs 
of winding numbers which count the times the trajectory 
winds around the torus in each direction, see Fig. ||) 

In view of the behavior shown in Fig. 3 and the re- 
sults of disorder diagrammatics, one may naively spec- 
ulate that the WL contribution to the structure factor- 
adds up in a similar way. Namely it consists of a series of 
singularities located at periods of the eight-shaped orbits 
illustrated in Fig. |l|c. One may also expect this contribu- 
tion to be positive, as in diffusive systems, since it should 
reflect an increase in the return probability compared to 
the classical value (i.e. the diagonal approximation). 

However, as we show here, this picture is inaccurate. 
Indeed, in the ballistic regime, some singularities do ap- 
pear at times which can be interpreted as periods of 
eight-shaped orbits (Fig. lc). But these contributions 
are rather weak. A large negative contribution comes 
from the original periodic orbits. It is superimposed on 
a smooth positive background which is not related to 
properties of the clean system. At certain cases the WL 
contribution to the structure factor can even become al- 
together negative. Thus, in ballistic systems, it does not 
have, necessarily, a definite sign. 

To make the paper self-contained, we organized it as 
follows: In the next section we prepare the mathematical 
background for our derivation by reviewing the standard 
results of disorder diagrammatics in the diffusive limit. 
This way we set the basis for extending the diagrammatic 
approach to the ballistic regime. In section III we 



FIG. 4. Periodic orbits of a particle moving on a tours are 
defined by pairs of integer numbers (n x ,n y ). These "wind- 
ing numbers" count the number of times the trajectory winds 
around the torus in the x and in the y directions, respectively. 
Some particular examples are: (a) the orbit (0,1) with length 
L, (b) the orbit (1,1) of length 2^L, and (c) the orbit (1,2) of 
length 

derive our central formulae for the WL contribution to 
R(uj), and the structure factor, S(t). In section IV, we 
analyze these results and derive an asymptotic expres- 
sion for S(t). Finally, we summarize and present our 
conclusions in Section V. 



II. BACKGROUND 



The purpose of this section is to lay the technical back- 
ground, and set the nomenclature for the analysis which 
will be carried out in the forthcoming sections. We shall 
review the main ideas-pf disorder diagrammatic technique 
for diffusive systemso, present the basic building blocks, 
discuss the approximations involved, and the limits of 
applicability. Finally, we summarize the results for the 
WL contribution to R(lo) within RMT framework, and 
for diffusive systems. These results will form a reference 
point for the analysis of i?(w) in the ballistic limit, which 
will be carried out in the next section. 

The disorder diagrammatic approach for Hamiltoni- 
ans of the type (ji|) is an efficient way of constructing 
the perturbation expansion, in the weak potential V(r), 
for quantities averaged over the disorder configurations. 
Examples of such quantities are n-point spectral corre- 
lation functions, the magnetic susceptibility, and various 
properties of the conductance of disordered metals. 

This diagrammatic approach is a semiclassical approx- 
imation in which the ratio of the particle wave length, 
Ai?, to the elastic mean free path, I, is assumed to be 
small. Therefore, it takes the formal form of an asymp- 
totic series in powers of I/UfI, where kp = 2n/Xp is the 
Fermi wavenumber. Yet, usually there will be also non- 
perturbative contributions, which are important when 
trying to resolve features on the scale of the mean level 
spacing, A, or over time scales longer than the Heisenberg 
time tn = 2irfr/A. Therefore, the applicability range of 
disorder diagrammatic is also limited to times smaller 
than the Heisenberg time, and energies larger than the 
mean level spacing. 

As a first example, consider the average of the retarded 
Green function: 
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FIG. 5. Example of diagrams contributing to the average 
Green function, Gf (k) (represented by the bold line). Thin 
lines represent the free Green function (i.e. in the absence of 
disorder), and dashed lines represents impurity scatterers. 



Gf (k) = 



IT) 



h 2 k 2 

2m 



V(r) 



Here (• • •) denotes an averaging over the configurations 
of the disordered potential, rj is an infinitesimal positive 
number, and p is the particle momentum. Expanding 
the Green function in powers of V(r), and changing rep- 
resentation to momentum space yields 

Gf (k) = G R (e, k) + G R (e, k) (V ) G R (e, k) + 



k' 



Gf (e, k)G«(e, k')G*(e, k) (F k _ k ^ k ,_ k ) 



where Vq = J <i 2 reR q ' r V(r) is the Fourier transform of 

the potential (|), and G R (e,k) = l/(e + ir)- is the 
free Green function. Terms containing an odd number of 
V-s vanish upon averaging, while those having an even 
number are calculated by Wick's theorem (since the po- 
tential V is Gaussian). Thus the average is equal to the 
product of averages of all possible pairs, such as (V^V_ q ). 
The various terms of this expansion can be represented 
diagrammatically as shown in Fig. |^. 

A partial summation of the infinite series of the dia- 
grams in Fig. ||, is achieved using Dyson's equation, and 
summation over the irreducible diagrams (those which 
cannot be separated into two disconnected diagrams by 
cutting one internal propagator line, e.g. (b) (d) and (e) 
m Fig. |). Thus the averaged Green function satisfies the 
relation 



Gf (k) = Gf (e, k) + Gf (k)SGf (e, k) 



(8) 



where £ is the self energy given by the sum over all 
irreducible diagrams, see Fig. [| To the leading order 
in l/kpl, £ is the contribution of the first diagram in 
Fig. |b. Thus 



S~^<V q F_ q )G (e,k + q) 
q 

hA 
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where P.V. denotes the principle value of the integral. 
The real part of £ can be absorbed into the definition 
of the reference energy e, thus the solution of Dyson's 
equation (||) yields 



FIG. 6. (a) The diagrammatic representation of Dyson's 
equation for the average Green function (^|). The bold and 
thin lines represents the dressed and the bare Green func- 
tions, respectively, (b) The self energy, E, given as a sum of 
irreducible diagrams. 



Gf(k) 



-c(k)- 
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where e(k) = (hk) 2 /2m. Similarly, the average of the 
advanced Green function is given by G^(k) = (Gf (k))*. 

Consider, next, the probability of a particle to arrive to 
r' in time t, when its initial state, |r; ep) , is a wave packet 
localized near r. We assume that this wave packet is com- 
posed of eigenstates centered at the Fermi energy, tp , and 
ranging over an energy band of width Ti/t. In the semi- 
classical limit Tl/t <C €f, these conditions imply that the 
particle velocity, vp , is well defined, and the wave packet 
width is of order of the elastic mean free path, I = vpT. 
The probability density for finding the particle at point 
r' after time t is given by P(r',r;i) = L 2 \U(r', r; t)\ 2 
where U(r',r;t) = (r'|e~^" fft |r; e^) is the propagator of 
the system. Using the convolution theorem, one obtains 



P(r',r;t) = ft due 



deV(r' , r; uj) 



(9) 



where 

V(r', r; u) = L 2 {G R (r', r; e F + hcj)G A (r, r'; e F )) , (10) 

and G R (r' , r; e) and G A (r' , r; e) are the exact Green func- 
tions of the system for particular realization of the dis- 
ordered potential. Notice that under our assumptions, 
2?(r', r; w) weakly depends on e, therefore the integration 
over e results in a factor of H/t. 

The diagrammatic expansion of T>(r',r;uj) proceeds 
along the same lines described above. It is convenient 
to perform the calculation in Fourier space, i.e. for 
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1 



dre 4q - r P(r',r' + r;w). (11) 



The leading contribution to 2?(q;cj), known as the dif- 
fuson, is given by the set of diagrams shown in Fig. 0a. 
The Dyson equation summing this set of diagrams yields 

2?(q,w) 



2wut 1 - n(w, q) ' 



where 



FIG. 7. The diagrams of the diffuson (a) and the Cooperon 
(b), and their interpretations as the contribution of pairs of 
classical orbits associated with the retarded and the advanced 
Green functions: (c) An orbit with itself (diffuson) and (d) 
an orbit with its time reversed counterpart (Cooperon). 
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1 1 



2tt £ 



fioj— hvp qcos6+Y- £- 
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(12) 



To obtain the second line of the above formula, we have 
expanded (k + q) 2 to linear order in q, and approximated 
2k ■ q by 2kpq cos 8, where Uf is Fermi wave number, and 
9 is the angle between the vectors k and q. This approx- 
imation is valid when gC/ip. 

In the diffusive limit, additional approximations can 
be made. Namely, one may use the small parameters 



ql <C 1 and lot <C 1, 



(13) 



to expand H(lo, q) in lot, 
form il(w, q) ~ 1 + iLOT — 
diffusion constant, thus 



and ql. The result takes the 
Dq 2 r, where D = 1 2 /2t is the 



(b) 




FIG. 8. The Hikami box associated with the interaction be- 
tween diffusons and Cooperons: (a) Its pictorial view in terms 
of "classical" trajectories, (b) Its diagrammatic expansion. 

be interpreted as the square root of the classical proba- 
bility to arrive to r', after time t, starting from r. Sub- 
stituting ( |l5| ) in (|l"rj|), yields Z>(r',r;w) as a double sum 
over the classical trajectories from r to r'. Approximat- 
ing the average of this double sum by the diagonal part, 
and substituting the result in (^) we obtain 



X>(q,w) 



lltVT —ilQT + Dq 2 T 



(14) 



This formula shows that the diffuson is the kernel of the 
diffusion equation: dn/dt = DW 2 n, where n(r) is the 
density of particles in real space. The diffuson is, there- 
fore, the classical mode of a disordered system in the limit 
of long time (lot <C 1) and large spatial scale ( ql <C 1). 

It isl. instructive to relate the diffuson to classical 
orbitscj. For this purpose we turn to calculate D(t,lo) 
using the van-Vleck approximation for the Green func- 
tions. A comment is now in order. The use of the van- 
Vleck propagator for a system with a white noise poten- 
tial is unjustified, since the scattering is not semiclassical. 
Therefore, here, we assume the disorder potential to be 
in the form of randomly located hard scatterers of size 
larger than the particle wave length. This potential is 
semiclassical, and produces diffusion on large scales of 
time and space. 

The van-Vleck's formula for the Green function, 
qR A^j-i r . a j s expressed as a sum over the classical 
trajectories 2 -!! from r to r' with energy e: 



G A (r, r'; e) ~ ^ B;e-^'^K 



(15) 



Here s /i (r',r; e) is the classical action of the /i-th trajec- 
tory, while is the corresponding amplitude which can 



where is the time which takes for the particle to travel 
from r to r' along the ^-th trajectory. Using classical sum 
rules, one can sum over the classical trajectories!?!! The 
result for diffusive systems is that of the diagrammatic 
calculation. This implies that set of diagrams associated 
with the diffuson is equivalent to the diagonal approxi- 
mation of pairs of orbits as shown in Fig. [7p. 

In systems with time reversal symmetry there is an ad- 
ditional classical mode called Cooperon. It comes from 
the infinite sum over the maximally crossed diagrams 
shown in Fig. ^b. These diagrams are obtained by revers- 
ing the direction of the momentum in one of the Green 
function lines. The classical picture of the Cooperon is, 
therefore, that of an orbit paired with its time reversed 
counterpart as shown in Fig. |?]d. It can be easily checked 
that the Cooperon has precisely the same analytical form 
of the diffuson. 

The issue of WL, in the language of diagrammatics, 
is the interaction between diffuson and Coopron modes. 
Pictorially, this interaction is the switching between the 
directions of the momenta of two trajectories, as shown 
in Figs, lc and 8a. The diagrammatic entity accounting 
for this switching is the Hikami boxE3, see Fig. ||b. It is 
a function of the incoming and outgoing momenta and 
frequencies of the diffuson and the Cooperon. For the 
particular choice of momenta and frequencies shown in 
Fig. ||b one has /i(q, q', lo) = h + h- + h + , where 
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h ° = E G *W k )G?+^( k -q-qO<^(k-q)Gf (k-q') 5 

k 

h- = ^—hi(q,-q';e,e + fi^)hi{-q,q';e,e + nw) (16) 
Zttt 2 2 

h+ = |— h\ (q, q'; e + e)/it (-q', -q; e + Tiw, e), 

2-7TT 2 2 

while 

hi(q, q'; 6, e') = £ G^(k)G*(k + q')G$(k + q). (17) 

k 

The calculation of the above diagrams in the diffusive 
limit ( |l3|) (the corresponding integrals are provided in 
Appendix A), gives 

47TT 4 

hM,u) = -^-[D{q 2 + q' 2 ) - ilo]. 

Having the basic ingredients of the disorder diagram- 
matics, we turn now to calculate the two-level corre- 
lation function defined by Eq. (|l|). Using the relation 
p{e) = Im{TrG fl (e)} /tt, we have 

R(u) = ^Re[(TrG fl (e + w)TrG A (e)) 

- {TrG R {e + Lo)){TrG A {e))] . 

This formula can be used as a starting point for diagram- 
matic expansion. However, it produces a large number of 
diagrams. A convenient way of reducing this number is 
to express R(u>) in terms of a generating function which 
has a simpler diagrammatic expansion. This generating 
function, r-F(u>), has been found by Smith, Lerner and 
AltshulerEJ. It satisfies the relation: 



R(u) 



If 



2^h 2 du- 



(18) 



and has the form of a free energy. The diagrammatic ex- 
pansion of F(to) can be loosely pictured as an expansion 
in the number of diffusons and Cooperons loops: 



F(w) = Fi(w) + F 2 {uj) + F s (w) + 



(19) 



Thus, the leading term, Fi(u>), is the contribution of 
the one loop diagram (see Fig. ||), i^C^O is the two-loop 
contribution (plus two additional terms whose role is to 
remove the ultraviolet divergence in the fpst diagram), 
F3(u) comes from three-loop diagrams, etco. In the pe- 
riodic orbit picture, Fi(u>) is the contribution of orbits 
shown in Fig. la+b, while F^w) is, in essence, the con- 
tribution of the eight-shape orbits illustrated in Fig. lc. 

The small parameter of the loop expansion ( |l9] ) is l/g, 
where g is the dimensionless conductance of the sys- 
tem, g oc tji/tc is the ratio of the Heisenberg time, 
tii = 2%%/ A, to the classical relaxation time of parti- 
cles in the system, t c . In diffusive systems, t c = L 2 /D 
(known as the Thouless timeEJ) is the time which takes 
for a classical particle to diffuse across the system. 







F 2 = 



FIG. 9. Diagrams of the free energy: Fi(u) is the leading 
contribution associated with the diagonal approximation of 
periodic orbit theory. Fi(u>) is the WL contribution to the 
free energy associated with the eight-shaped orbits of Fig. lc. 
Dashed impurity lines represent large momentum transfer, 
k > l/l, dotted lines represent small momentum transfer, 
k < l/l. 

The form of the free energy (|l9| ) together with Eq. ( |lg| ) 
induces a similar expansion for the two-level correlator: 

R(u) = Ri(w) + R 2 (oj) + Ra(u) + • • • , 



where 



Rj(w) 



A 2 d 2 
2^h 2 dcu 2 



ReF.H, j = 1,2,3 ■ 



(20) 



Thus Ri(ui) is the result of diagonal approximation, 
i?2(w) is the WL contribution, and additional terms give 
higher WL corrections. 

The leading contribution to the two-level correlation 
function, i?i(w), haa-been discussed extensively by Alt- 
shuler and ShklovskiO. It is straightforwardly calculated 
using (|20|). Taking into account the l/n symmetry factor 
of the n ladder diagram defining Fi(u>) (see Fig. 9) we 
obtain Fi(lu) = — ln(_Dg~T — itor), where the diffu- 
sive approximation (O) has been assumed. Notice that 
although this sum does not converge, its second deriva- 
tive with respect to u> does. Moreover, one can check 
that, in two dimensions, Ri(u>) = for u) > 0. Thais the 
leading term in this case is the WL contributions. 

In this paper we focus our attention on the WL contri- 
bution to R(uj) of two dimensional ballistic systems. As 
a reference point, however, it will be instructive to re- 
view results of RMT, and disorder diagrammatics in the 
diffusive limit. In both cases, our starting point is the 
diffusive form of the WL contribution to the free energy 
(obtained from the diagrams shown in Fig. 9): 



hir ^ (Dq 2 



iu))(Dq' 2 — iu) 



(21) 



The RMT result corresponds to the zero mode con- 
tribution (q = q' = 0) in the above sum, namely 
i*2 (u) ) ~ —iA/h'KUj. It is purely imaginary, therefore, 
Eq. (EOh implies that R 2 (u>) = in the RMT limit. Since 
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RMT accounts for the universal behavior of chaotic sys- 
tem, we conclude that R 2 (lo) is a purely nonuniversal 
quantity. 

Turning to the diffusive limit, we first note that the 
sum in (Ell) diverges logarithmically, even after differen- 
tiating twice with respect to u. Thus one has to intro- 
duce an upper cutoff on the momentum, which is usually 
taken to be 1//, where I is the elastic mean free path. As 
will be shown in the next chapter this artificial cutoff can 
be avoided if the approximations ( |l3| ) are not used in the 
calculations. 

To evaluate F 2 in the regime l/t c <C 1/r, one 

can also use dimensional regularizationli.3: Replacing the 
sums over q and q' by integrals, and evaluating them in 
d = 2 + r\ dimensions yields 



d d q 



1 



(2n) d Dq 2 



(22) 



Changing variables from q to (iD /tof-^l and using the 
formula (see Appendix A) 



J l + q 2 



one arrives at 



F 2 (u) 



/wr(47rL>) 2 \AttD ) V 2 ' 



R 2 {u)) is now obtained by taking the second derivative 
with respect to as follows from Eq. (ppf). Thus using 
T(l + 7}/2)T(-T]/2) = tt/ sin(-?77r/2) we have 



A7r(l + r))r)Re[(-iui) 



77-11 



2hg 2 (4irD)v sin 2 (-7 ?7 r/2)r 2 (l + ? ? /2) ' 



where g = An 2 T%D / L 2 A is the dimensionless conductance 
of the system. Finally we let r\ — > 0, and obtain 



R(u) ~ R 2 (w) 



A 



g 2 hu) 



1 1 

— < uj < -. 

tr. T 



Note that the domain of validity of the above formula 
vanishes in the ballistic limit since the classical relax- 
ation time, t c , is smaller or equal to the scattering time, 

T. 



to the leading order result in the small parameter l/L 
(since ql <C 1, q is of order l/L, and I -C L). In the 
ballistic regime this approximation cannot be used, and 
one has to evaluate integrals, such as ([l2|) and (16), to 
all orders in l/L. Moreover, diagrams having a small 
number of impurity lines form the dominant contribu- 
tion (unlike in the diffusive regime), therefore, possible 
cancellations among diagrams as well as double counting 
should be examined carefully. The outcome of this exam- 
ination is that diffusons and Cooperons contributing to 
F 2 (to) should start from two impurity lines. Apart form 
this point, the WL contribution is given by the same di- 
agrams shown in Fig. 9, but evaluated to all orders in 
l/L. 

We begin by deriving an expression for the diffuson 
(starting from two impurity lines) in the ballistic regime. 
Dyson's equation, in this case, yields 



n(^,q) 



2-kvt 1 - Jl(<j, q) ' 



X>(q,w) 



where n(w, q) is the integral given by Eq. (|12[). To avoid 
the expansion in ql and cut, here we first integrate over 
£ (by closing the contour in the complex plain), and 
then integrate over the angle, 9, exactly. The result is 
II(w,q) = 1/Q u (q), where 



Q u (q) - ^{l-iLOT) 2 + {lq) 2 . 
Thus the generalized formula for the diffuson is 



h 



2ttvt Qu{q) - 1 



(23) 



A similar calculation for the Cooperon produces the same 
analytical expression. 

The above formula holds both in the ballistic as well 
as the diffusive regime dl3| ) . It can be easily checked that 
an expansion of the denominator of Eq. ( |23|) in lot and 
ql, yields the result of the diffusive limit, (|l4|). 

The calculation of the Hikami box (Fig. 8b, Eq. |l^), in 
the ballistic limit, follows along the same lines. Namely, 
one first integrates over the modulus of k, and then the 
remaining angular integration is performed exactly. For 
example, after integration over the modulus of k, Eq. (|lT 
reduces to an integral of the from 



III. WEAK LOCALIZATION IN THE 
NONDIFFUSIVE REGIME 



1 f 271 

2ttJ (l+xi 



dO 



cos<9)(l + a;2 cos(# — tp)) 



(24) 



In this section we calculate the WL contribution to the 
two-point correlator in the ballistic regime. By ballistic 
we refer to the situation in which the elastic mean free 
path, I, is of order of the size of the system L. To under- 
stand what kind of changes are needed in order to extend 
the diagrammatic calculation into the ballistic regime, 
recall that the diffusive approximation (n~3) corresponds 



The result of the integration over the angle 9 (see Ap- 
pendix A) is 



I(x 1 ,x 2 ,(f) = [ — + — 



1 



Vi 2/2/ I + yiy 2 - x^ 2 cos(ip) 



(25) 



where 
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FIG. 10. The results for the WL contribution to the struc- 
ture factor at various values of the ratio of the elastic mean 
free path, I, to the size of the system L. The WL effect, in 
this system, becomes stronger in two limits: (a) the diffu- 
sive regime, I <C L, where the particle approaches localization 
in real space, and (b) the ballistic limit, I 3> L, where the 
particle becomes localized in momentum space. 



FIG. 11. The weak localization contribution to the struc- 
ture factor, in the ballistic regime, exhibits a singular behav- 
ior. The singularities are located at times which are linear 
combinations of periods of two orbits of the clean system (see 
Fig. 4). Here t(n x ,n y ) = (n% + riy) 1 ^ 2 L/vf denote the period 
of the periodic orbit defined by the pair of winding numbers 
(n x ,n y ). 



1,2. 



With the help of this function, the various terms of the 
Hikami box (see Eq. ([!(]) and Fig. 8b) take the form: 



2ttt 3 



ft d (l-iwr) 4 A 



ilqi 



±ilq 2 



1 — iu)T 1 — ilot' 



¥>12 



where ipi 2 is the angle between qi and q 2 , and 



47tt 3 (1 — iwr) 



2/7 \ 2 
/ I-Q2 



(Jqi 
\V2fJ 'VV2/2 



where 



fi = V(l - iur) 2 + (*072 i = l,2. 
Collecting the diagrams of ^2(0;) (Fig. 9) we obtain 



2 A 2 



471-2x2 



1,2 



where we use the notation 

T>i — 2itvTT){<\i,u)) /h, 

and the sum is over the vectors qi and <\ 2 . The pe- 
riodic boundary conditions in our system imply that 
q; = 2irm.i/L where is an integer vector of two com- 
ponents. 

The above formula is our central result. Performing the 
sum over momenta and substitution it in ( pt]| ) gives the 
exact WL contribution to the two-level correlation func- 
tion in the semiclassical limit. The applicability range 



of our result goes beyond the diffusive limit, I -C L, and 
includes the ballistic regime, I ~ L, as well. In contrast 
with the formula in the diffusive limit (Eq. ^l|), here the 
momenta sum converges, and there is no need to intro- 
duce an arbitrary cutoff or regularization. The results 
which will be shown below were obtained by performing 
the momenta sum numerically with cutoff chosen such 
that contribution of additional terms is of order of the 
numerical error. 

In presenting our results it will be convenient to em- 
ploy the spectral structure factor defined in Eq. (0). We 
denote by S 2 (t) the corresponding WL contribution, 



A 



dioR 2 (u>)e 



(26) 



and rescale its magnitude by a factor of 2ir 3 g 2 , where 
g cx tn/tc is the dimensionless conductance of the sys- 
tem. Note that in the ballistic regime, the relaxation 
time, t c , is no longer the diffusion time. It is approxi- 
mately the traversal time across the system, t c = L/vp 
where vf is the velocity of the particle, and L is the size 
of the system. Therefore, from now on we define g to be 



9 = 



Tivf 
AL' 



(27) 



In Fig. [h] we plot S 2 (t), for various values of the ratio 
between the elastic mean free path and the size of the 
system. These values range from diffusive (l/L = 0.01) 
to ballistic (l/L — 1.5) dynamics. Several features of 
S 2 (t) are evident: First, the WL contribution appears 
only within a finite interval of time. It vanishes both at 
t = and when t — > 00. Second, in both limits, I <C L 
and I 3> L, the WL contribution becomes strong. Third, 
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A //L=0.5 





















FIG. 12. The convergence behavior of the momenta sum of 
the structure factor in the ballistic regime. The smooth part 
(dash-dotted) is determined by the zero mode, qi = q2 = 0. 
The dashed line is the result of a sum over momenta within 
radius of 2.837T. Higher q terms build up singularities along 
the time axis as demonstrated by the solid line where the 
momenta sum extends to radius of SOn. 

in the ballistic regime, I ~ L, S 2 (t) exhibits a distinctive 
singular behavior consisting of a series of dips. These 
dips are located at times which are combinations of peri- 
ods of the periodic orbits of the clean system. In Fig. |ll| 
we depict 5*2(0 f° r I = L/2 and indicate the singularities 
with the corresponding combinations of periodic orbits. 



IV. ANALYSIS 

In analyzing the above results it is instructive to study, 
first, the convergence behavior of the momenta sum of 
the WL contribution. In Fig. [l2] we depict the results 
for 5*2 (t) calculated in the following approximations: The 
dash-dotted line is the contribution coming from the zero 
mode, qi = q2 = 0. Clearly this mode dictates the gross 
behavior of S 2 (t). In particular, it determines the inter- 
val of time where WL is significant. The dashed line is 
the result obtained by taking into account the next lowest 
momentum modes, i.e. summing over qi and q2 within 
the radius qi,q2 < \/8tt. In this approximation some ad- 
ditional features of S 2 (t) are resolved. The solid line is 
the result of the full momenta sum. Thus the singular 
behavior of the structure factor comes from the tail of 
the sum. 

To obtain a simple analytic characterization of the WL 
contribution to the structure factor, we proceed in the 
following way. First, we derive a formula for the smooth 
part of 5*2 (t) given by the contribution of the zeroth mode 
qi = q2 =0. This formula will give us the main param- 
eters characterizing the WL contribution in the ballis- 
tic limit. Then, we evaluate the momenta sum in the 
asymptotic limit of large u>. The result of this calcula- 
tion provide the local behavior of S 2 (t) in the vicinity of 



the singularities. 

To calculate the smooth part of 5*2 (t), denoted here- 
inafter by S2 (t), we start by evaluating the zero mode 
contribution to #2(0;). A straightforward calculation of 
the term qi = q 2 = yields 



R 2 (u>; q = contrib.) = 



4r 3 A 3 [5_54( W t) 2 + 21(wt) 4 ] 



h 3 TT 3 [l + (ujT) 2 } e 
Taking, now, the Fourier transform we obtain: 
I 2 



s 2 {t) 



\2it 2 g 2 L 2 



rif + mr + 6t + 6), 



(28) 



where t — t/r. We remark, here, that the above for- 
mula applies only in the ballistic regime, where the ze- 
roth mode is dominant. In the diffusive systems, the zero 
mode contribution is negligible compared to that coming 
from the momenta sum in regime 1/L < q <C l/l (which 
is absent in the ballistic case). 

Formula (pit) allows one to characterize the major fea- 
tures of the WL contribution to the structure factor: The 
time t* where S 2 (t) is maximal; its value at this point, 
S 2 ; and the width of the time interval where the WL 
effects are appreciable, W*. The results are: 

t* = 4.24 r, 
W* = 2.46 r, 



and 



S.; = 1 ).:-!">:-! ( — 
Lg 



0.:-;.-:-; I 1— 



where the interval width is defined by (W*) 2 = 
J dtS 2 (t)(t - t*) 2 / J dtS 2 (t). Thus WL effects, in the 
ballistic limit, are pronounced within a time interval of 
width 2.46t centered at t = 4.24r, and the typical value 
of the WL contribution is proportional to r 2 . 

Note that these results are independent of the size of 
the system. Therefore the gross behavior of S 2 (t) is not 
influenced by the periodic orbits of the clean system. It 
is natural to ask what is the role the classical orbits of 
the system? As we show below, these orbits lead to the 
singular features decorating the smooth part of the struc- 
ture factor S 2 (t) as demonstrated in Fig. 12. 

In analyzing this singular behavior, we first notice that 
its main part comes from large uj or equivalently large val- 
ues of the momenta qi and q2. Therefore, to calculate 
this contribution it is sufficient to approximate the dis- 
crete angular sum of F 2 (lo) (over the phase between the 
vectors qi and q2) by an integral. The small parameter 
of this approximation is 1/wr. From (|2g) one finds that 
the angular average of the WL contribution to the free 
energy, denoted by F 2 {uj), is: 



F 2 {io) 



rA 



A + B + C 



^ QMQi 

qi,q2 v 



Q2XQ1-IXQ2-I)' 
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where 



A = 



1 — iujt ' 



D 



and 



C 



-B(Q 1 + Q 2 ), 



At asymptotically large values of lot the leading contri- 
bution comes from C . This is evident once noticing that 
when cut — > oo, Qi — ► lut, and therefore A, B = 0(— ) 
whereas C = 0(1). 

Next we apply the Poisson summation formula to con- 
vert the sum over qi and c\ 2 into an integral. The free 
energy is then expressed as 



(m,n) 



(29) 



(n + m)L/vp, and 0(x) is the step function. The am- 
plitude of each contribution, B nm , depends on the val- 
ues n and m. For cases where either n or m vanish, 
Bom = B n fl ~ e 3 / 2 /12\/27r 4 , while if n and to are large, 
B nm ~ e 3 /144?i 2 7r 4 . 

Thus, S 2 (t) is composed of a smooth contribution 
(E8fi and a sequence of singular functions, of the form 



-re 



-(t-tm»)/r 



0(t — t rnn ), where i m „ is the period of a 
composite orbit, i.e. the sum of periods of two periodic 
orbits of the clean system. Each singular contribution is 
negative, and its magnitude at time t mn is proportional 
to t 2 nn e~* m "/ T /Z 4 . The contribution associated with sin- 
gle orbits (i.e. when either norm vanish) is considerably 
larger than that of composite orbit (in which both n and 
to differ from zero). In any case, the singular contribu- 
tion decreases exponentially in time , and as a power law 
in 1/7 (when L < I). This behavior is indeed observed in 
Figs. 10, 11, and 12. 



where m and n are integer vectors. As will be shown 
below, these integer vectors are associated with winding 
numbers of the periodic orbits of the clean system. Each 
term in formula (p9|) is of the form 



where 



K(m) 



tAL 4 



-(1 - iujT)K(m)K(n), 



dq qJo(mqL) 
o 27 Ql(q)(QUq) - 1) 



(30) 



(31) 



Here Jq{x) is the Bessel function of zero order, and 
to = |m| is the magnitude of the vector m. For to = 
this integral yields 



K(0) 



1 



271-/2(1 - iujT 



+ 0(l/c 2 ) 



For m 7^ the integral ([31]) can be calculated using the 
steepest descents method (see Appendix B), and in the 
large to limit it gives 



K(m) 



=3/2 



r) 



72ttZ 2 



1 — IUJT 



(32) 



The above results imply the following form of the struc- 
ture factor 

s 2 (t)~s 2 (t) + Y,st m \t), 



where 



Si n ' m \t) = -B n 



~ ) B[/ - /,„„] Pr 



is the contribution associated with orbits characterized 
by the winding vectors n and m. Here t — t/r, t nm = 



V. SUMMARY AND CONCLUDING REMARKS 

In this paper we have calculated the WL contributions 
to the two-level correlation function and its Fourier trans- 
form, the structure factor. These are the leading quan- 
tum interference effects which affect the spectral statis- 
tical properties of the system defined in (Q) . 

Our theory generalize previous calculations of the 
WL coatribution to the spectral statistics of diffusive 
systemscJ by extending them into the ballistic regime 
where the elastic mean free path, I, is of order of the size 
of the system, L. Here the disorder is weak enough to 
leave traces of the dynamics of the underlying clean sys- 
tem, which appear as singularities in the structure factor 
(Figs.0, 0, and©. 

Our study has been focused on spectral rather than dy- 
namical characteristics to avoid the problem of specifying 
the manifold on which WL takes place. Indeed Fig. 
demonstrates that the WL contribution is pronounced in 
two limits. Panel (a) of Fig. 10 is a representative ex- 
ample of the results deep in the diffusive regime I <C L, 
while panel (b) shows the typical behavior in the ballis- 
tic limit, / = 1.5L. In both cases, the system approaches 
the strong localization limit, but the localization is of dif- 
ferent nature. In the diffusive case, it is localization in 
real spacea, whereas in the ballistic case the localization 
is on a quasi-one-dimensional annulus in the momentum 
space. (This is evident once noticing that on clean torus, 
eigenstates are plain waves and therefore the particle is 
localized in momentum space.) In the latter case, it is 
suggestive that the effective dynamics is associated with 
Levy flightsEj rather than diffusion, since the disorder 
couples, predominantly, momentum states with degener- 
ate eigenvalues, .which may lie far away along the mo- 
mentum annulusEZI. 

A simple semiclassical interpretation of our results, 
within periodic orbit theory, is not straightforward. The 
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results, clearly, cannot be obtained from a diagonal ap- 
proximation in which higher order fi corrections are 
added to Gutzwiller's trace formula (e.g. diffracting or- 
bits, creeping orbits, etc. ). One can easily verified that 
such approximation yields only a positive contribution, in 
contrast with our results. This type of correction might 
explain the smooth positive part of the WL contribution, 
§2(t). However, a correct analysis within periodic orbit 
theory must go beyond the diagonal approximation, and 
take into account pairing of orbits which are not related 
by symmetry, but have actions exponentially close one to 
the other (up to a constant phase, 7r, which is needed in 
order to explain the negative contribution of the periodic 
orbits). The fact that the WL contribution may become 
negative at certain times implies that the system exhibits 
anti-weak-localization at certain regions in phase space. 
This may be related to anti-scarring effects observed in 
wave functions of chaotic billiardsea. 

Nevertheless, our work still elucidates several features 
of the leading WL effects in ballistic chaotic systems. 
First it shows that it appears within a finite interval of 
time; Second, it has a singular behavior associated with 
periodic orbits and linear combinations of periodic or- 
bits; Third it can have different signs at different points 
in phase space. 

These results have important consequences: First they 
show that the dominant contribution to the WL, in the 
ballistic regime, does not come from the eight-shaped 
orbits (Fig. lc), as suggested by the diagrammatic pic- 
ture. The main contribution comes from diffracting 
orbits (which are not related to the classical periodic 
orbits of the system), as well as from the original pe- 
riodic orbits of the system. Moreover, the zero mode 
contribution, defining ^(t), plays a dominant role here, 
while according to the results of the ballistic er-modcl 
it should vanish (since the zero mode of the cr-model is 
identical to RMT). The apparent contradiction between 
our results and those of the ballistic cr-model is proba- 
bly due to the fact that the ballistic cr-model does not 
account correctly for the return probability. This is also 
manifested by the so called "repetition problem" which is 
a small mismatch, associated with repetitions of periodic 
orbits, between the exact asymptotic results of periodic 
orbit theory and those of the ballistic cr-model. Ideas 
associated- with memory effects in long range random 
potentials may be useful in resolving this problem. 
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Appendix A. Useful integrals 



In this appendix we calculate useful integrals fre- 
quently encountered when calculating diagrams which 
appear in this paper. The first type of such integrals 
appear when integrating products of retarded and ad- 
vanced Green functions over the energy. The integral is 
of the form: 



Y 

1 n,' 



dr\ 



1 



1 



n 



where n and m are non-zero integers. Applying the 
Cauchy theorem, and using the fact that the co- 
efficient a_i of a Laurent series, ^2i<i-i(z — zo)~ l 
of a function with an n-th order pole is a_i = 
[{z — Zo) n f(z)] , one immediately gets 



(n-l)! dz 



y 

1 71,' 



27t(to 



2) 



| -rn— n m+n— 1 



(n- l)!(m- 1)! 



The second type of integral appears when integrating 
over momenta in d dimensions, e.g. in the calculation 
of the WL contribution to the free energy ( p2| ) in the 
diffusive limit. This family of integrals is of the form 



L n,d 



d d q 



(l + q 2 ) 



2\n ' 



where n is an integer, d is a real number, and d d q = 
dVLq d ~ 1 denotes the measure in d dimensions. Since the 
integrand is independent of the angles, the angular inte- 
gral yields fdfl = dTT d / 2 /T(l + d/2). This formula should 
be understood as an analytic continuation of a function 
defined on an infinite set of integer values of d. Changing 
the integration variable to x = q 2 yields 



In 4 
■nid 



TT d ' 2 d 



2r(i + rf/2) 

n— 1 Qn— 1 



dx- 



r d/2-l 



(-1) 



2r(l + f) (n-l)W(3 n - 



(1 + x) n 
d£ I dxx^ 1 



Integrating over x, then changing the integration variable 
to y = and integrating over y we obtain 



7r d / 2 r(d/2)r(i-rf/2)(-i)"- 

T(n)T(l + d/2-n) 



As the last step, we use the relation T(x + 1) 
simplify the expression. The result is: 



xT(x) to 
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d/2 T(n-d/2) 

In ' d -* r(n) ■ 

We turn now to evaluate the integral which appears 
when calculating diagrams in the ballistic limit, namely 
the integral defined by Eq. (|24|). It will be calculated in 
the regime \x\\, \x^\ < 1 (corresponding to the diffusive 
limit) and then analytically continued to the full complex 
plain. It is natural to substitute z = e lS which immedi- 
ately transforms (|24|) into the contour integral 



I(xi,x 2 ,tp) 
-2i 



zdz 



7TXiX2 7j z |=i (z 2 + £z+l)(z 2 e- i 'f> + JLz+e i! f>) 
-2ie iv f zdz 



TTXIX2 J\z\=l ( z - z l)( z - Z 2){Z - Z 3 )(Z - Zi) 

where Zi are 



■23 




z 2 



Zi 



1 

Z\ 
Zs 



= — (33) 



Only the poles z\ and 2:3, which lie inside the unit cir- 
cle, contribute to the integral. Thus using the residue 
theorem we get 



I(xi,x 2 , if) 



X1X2 



(Zl - Z 2 ){Z1 - Z 3 )(Z1 - Zi) 

Z3 



(34) 



(^3 - Zl)(Z3 ~ z 2)( z 3 - Zi). 

Finally, we substitute ( |33|) in ([}4]) and arrived at (^ 
Appendix B. Asymptotics of K[m) 



In this appendix we evaluate the integral ( |3l"| ) in the 
asymptotic limit m 3> 1 and lot 3> 1. We begin by 
changing variables to r, — rrfj^ and taking the leading 



term in l/(i + lot) . The result is 



K(m) 



2wl 2 (i + lot) 



{rf — 1)2 



Being interested in the leading order expansion in 1 /lot, 
we further approximate the integral by substituting the 
asymptotic formula of the Bessel function: Jo(x) ~ 

y / ^"cos(a; — 7r/4) as x — » 00. Representing the cosine 
as a sum of two exponents we arrive at: 
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FIG. 13. The contour of integration for the asymptotic 
evaluation of K(m). It is composed of steepest descent paths 
(a) and (b), connected by arcs, (c) and (d), which give van- 
ishing contributions. The dashed lines represent cuts of the 
integrand in the complex plain (see Eq. Bq) 



K{m) 



+ LOT)° 



x > / — : — : — 3-e 

(rj 2 — 1) 2 



(35) 



The two terms in the above sum will be handled sepa- 
rately. Later, it will become clear, that the leading con- 
tribution to the integral comes from the plus-sign term. 
Therefore, for the time being, we ignore the term with 
the minus sign. Using the Cauchy theorem we can de- 
form the contour of the integrals such that its direction 
near the edge at 77 = is a steepest descent direction. 
The contour is further deformed to follow steepest de- 
scent curves as illustrated in Fig. [ll| Thus the contour 
consists of four segments: (a) from to ioo/(l + i/lot), 
(b) the part surrounding the cut, (c) an arc connecting 
these two segments at infinity, and (d) an arc connect- 
ing the end of the (b) path and the original end point at 
(lot — i)oo. 

It is straightforward to check that the contribution 
from part (a) is of order 1/(wt) 7 / 2 . This contribution 
will turn out to be negligible compared to the one coming 
form segment (b). It is also clear that the contributions 
of arcs (c) and (d) vanish, when their distance form the 
origin approaches infinity. Thus we focus our attention 
on segment (b). 

To evaluate the integral we exponentiate the pre- 
exponent factor in Eq. (|35|), and write the integral in 
the form 



K(m) 



+ lotV 



where 



A{r,) 



i— — (i + LOT)rj — — \n(r] — 1). 



The steepest decent contour is found in the usual way: 
First, one takes the derivative of A{rf) with respect to 77, 
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A \ri)=i— (i + ojt -3/2 -, 

i 77 — 1 

and find the saddle point 77*, satisfying A' (77*) = 0. The 
result is 



77* = u* + if* = 1 — 



3/ 



2mi(l — ilot) 



where u* and v* are real numbers. Second, one substi- 
tutes 77 = (u — u*) + i(v — v*), where u and v are real, 
and solve for the curve which satisfies the condition 

lm(A{rj)) = lm(A(rf)) 

The formula for this curve is 



v — uiu 2 

= tan — (uiu — v) . 

u -\- uiv — 4 3 



Rotating the axis as x 
obtains 



uiu and y = u + uiv, one 



^ 2 tan %x 



This exact form shows that the contour can be indeed 
deformed as shown in Fig. 13. It also provides the possi- 
bility of constructing the full asymptotic expansion of the 
integral. However, in view of the approximations which 
already have been made, we are interested only in the 
leading term. 

Thus, taking the quadratic approximation for the ac- 
tion: A(rj) ~ A(r)*) + \A"{rf) (77 - 77* ) 2 , with 

A"(v) = l, ' 



2 (77- l) 2 ' 

and evaluating the resulting Gaussian integral we arrive 
at ([32|). Notice that the result is proportional to \/lot. 
Thus, the edge contribution (which is of order l/(wr) 7/2 ) 
can be indeed neglected. This calculation also shows that 
A" (rf) oc to 2 , and therefore the small parameter of this 
saddle point approximation is 1/to. 



The contribution of the second term in the sum (35), 
i.e the one with the minus sign, is calculated following 
along the same lines. However, it turns out that in this 
case the deformed contour does not pass through a sad- 
dle point, and the only contribution comes from the edge 
at 7/ = 0. It is, again, of order 1/(wt) 7/2 , and therefore 
can be neglected. 
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